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We report a search for B — > J/ip D decays, based on a sample of 124 x 10 6 BB events collected 
with the BABAR detector at the PEP-II storage ring of the Stanford Linear Accelerator Center. 

No significant signal is found. We obtain upper limits on the branching fractions of 1.3 x 10 -5 
for B° -> J/ipD° and 1.2 x 10 -4 for B + -> J/ip D + at 90% confidence level. 

PACS numbers: 13.25.Hw, 14.40.Nd 



Measurements of the inclusive spectrum of charmo- 
nium mesons in B decays are in conflict with conventional 
expectations. The spectra of the momentum of the Jhb 
mesons in the TY45) rest frame observed by CLEO [lj 
and by BaBar |2j (Fig. [IJ, compared with calculations 
using non-relativistic QCD (NRQCD) show an excess 
at low momentum, corresponding to a branching fraction 
of approximately 6 x 10~ 4 . Various hypotheses have been 
proposed to explain this low-momentum excess. 

Brodsky and Navarra Q have suggested that the de- 
cay B —>■ J/ip Ap Q , with the possible formation of a A-p 
bound state, could explain the CLEO result. The kine- 
matic boundary of this structure corresponds to the case 
where the J/ip recoils nearly monoenergetically in the B 
rest frame against a 2GeV/c 2 particle. The A-p state 
could be observed near or just below threshold. BABAR 
has searched for these decays and obtained an upper limit 
of 2.6 x 1CT 5 at 90% confidence level (C.L.) Q, too small 
to support the mechanism proposed in Q. 

Decays to a J /ip meson and a hybrid meson, i.e. a 
bound state of two quarks and a gluon, have been pro- 
posed 0,0- I n this case the hybrid meson, possibly a 
(sdg) state, would need to have a mass of about 2 GeV/ c 2 . 
No experimental evidence has been found to support this 




0.5 1 1.5 2 



PCM ( GeV/c ) 

FIG. 1: Upsilon(4S) rest-frame momentum of J/ip mesons 
produced in B decays, after subtracting feed-down from 
Xci.a -» J/tpl and ip(2S) -> J/ipw-K (points) 2]. The 
histogram is the sum of the color-octet component from a 
NRQCD calculation Q (dashed line), which includes multi- 
body final states, and the color-singlet J/ipK^*> component 
(dotted line). The normalization of the curves has been con- 
strained to fit the data. 



mechanism. 

If B mesons were decaying to a narrow resonance and 
a J/ip meson, the J/ip meson would be monoenergetic in 
the B rest frame. Such peaks would appear smeared with 
an RMS of 0.12 GeV/c in Fig. ID due to the motion of the 
B in the T(4S) rest frame. 

The presence of bucc components (intrinsic charm) in 
the B-meson wave function has also been proposed. In 
that case the charmonium meson is obtained merely by 
dissociation when the b quark decays. Intrinsic charm 
was first introduced by Brodsky et al. to explain 
an unexpectedly large cross-section for charmed-particle 
production in hadron collisions. Using the estimated 
amount of intrinsic charm in the proton as an input, 
Chang and Hou predict B — ► J/ipD(w) decays with 
branching fractions of the order of 10 -4 [lOj. The dom- 
inant final state is expected to be B — ► J/ipDir, for 
which BABAR has reported an up per limit at 90% C.L. of 
5.2 x 10~ 5 for B+ -> J/ipD°Tr+ 11]. Four-body decays 
such as B — > J/ipDirir should be extremely suppressed 
by the small phase space available near the kinematical 
limit. The remaining untested final state is J/ipD. Cal- 
culations by Eilam et al. using perturbative QCD [l^ 
predict branching fractions (BF's) for B — > J/tpD decays 
on the order of 10 -8 - 10 -7 . The observation of a signal 
with a BF significantly larger would suggest the presence 
of intrinsic charm inside the B meson. 

In this Letter we report a search for decays B — > J/ipD, 
with D° decaying to K+tt^, D+ to K°n+ , K° to tt+tt", 
and J/ip — > £ + £~ , where £ is e or \i. 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II storage ring and comprise 
an integrated luminosity of 112 fb _1 taken at the T(4S) 
resonance. The BABAR detector is described in detail 
elsewhere jj^- A five-layer, double-sided silicon vertex 
tracker (SVT) surrounds the interaction point and pro- 
vides precise reconstruction of track angles and B-decay 
vertices. A 40-layer drift chamber (DCH) provides mea- 
surements of the transverse momenta of charged parti- 
cles. An internally reflecting ring-imaging Chcrcnkov 
detector (DIRC) is used for particle identification. A 
CsI(Tl) crystal electromagnetic calorimeter (EMC) de- 
tects photons and electrons. The calorimeter is sur- 
rounded by a solenoidal magnet providing a 1.5-T field. 
The flux return is instrumented with resistive plate cham- 
bers used for muon and ncutral-hadron identification. 
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We select multihadron events by demanding a mini- 
mum of three reconstructed charged tracks in the polar 
angle range 0.41 < 9^ < 2.54 rad. A charged track 
must be reconstructed in the DCH, and, except for the 
reconstruction of —>■ ir + ir~, it must originate at the 
nominal interaction point to within 1.5 cm in the plane 
transverse to the beam and to within 10 cm along the 
beam. Events are required to have an T(AS) production 
point within 0.5 cm of the average position of the inter- 
action point in the plane transverse to the beamline, and 
within 6 cm longitudinally. Neutral clusters are defined 
as electromagnetic depositions in the calorimeter in the 
polar angle range 0.410 < 6i a b < 2.409 rad that are not 
associated with charged tracks and that have an energy 
greater than 30 MeV and a shower shape consistent with 
a photon interaction. We require the total energy for 
charged tracks and photon candidates in the fiducial re- 
gion to be greater than 4.5 GeV. To reduce continuum 
e + e~ — * qq background, we require the ratio of second- 
to-zeroth Fox- Wolfram moments i?2 0] of the event, cal- 
culated with both charged tracks and neutral clusters, to 
be less than 0.5. Charged tracks are required to be in re- 
gions of polar angle for which the particle identification 
(PID) efficiency is well measured. For electrons, muons, 
and kaons the acceptable ranges are 0.40 to 2.40, 0.30 to 
2.70, and 0.45 to 2.50 rad, respectively. 

We further select signal events as described in the fol- 
lowing. Event selection is optimized by maximizing the 
sensitivity s = e/ (ci/2+a/ZYb), where a = 3 is the number 
of standard deviations of significance desired 0. The 
maximum of this ratio is independent of the unknown 
signal branching fraction. The signal efficiency e after all 
selection requirements is estimated from simulated Monte 
Carlo (MC) samples. The number of background events 
Nb, scaled to the integrated luminosity of the data, is 
estimated using inclusive T(AS) — > BB and e + e~ — > qq 
MC samples. 

We reconstruct J/tp candidates from a pair of oppo- 
sitely charged lepton candidates that form a good ver- 
tex. Muon (electron) candidates are identified with a 
neural-network (cut-based) selector. For J/tp — > e + e~ 
decays, electron candidates are provisionally combined 
with nearby photon candidates in order to recover some 
of the energy lost through bremsstrahlung. These 
bremsstrahlung-photon candidates are characterized by 
a deposit of more that 30 MeV in the electromagnetic 
calorimeter and a polar angle within 35 mrad of the elec- 
tron direction, as well as an azimuthal angle either within 
50 mrad of the electron direction, or between the electron 
direction at the origin and the azimuth of the impact 
point in the EMC. The lepton-pair invariant mass must 
be in the range [3.00, 3.14] GeV/c 2 for both lepton fla- 
vors. 

We form K® candidates from oppositely charged tracks 
originating from a common vertex and having an invari- 
ant mass in the range [487, 510] MeV/c 2 . The K° flight 



length must be greater than 1 mm, and its direction in 
the plane perpendicular to the beam line must be within 
0.2 rad of the K® momentum vector. All charged tracks 
are taken as pion candidates, and kaon candidates are 
identified with a likelihood selector based on Cherenkov- 
angle measurements from the DIRC and specific ioniza- 
tion in the SVT and in the DCH. Candidates for D 
mesons are formed from Kn combinations; a requirement 
on the Kn invariant mass ttik-k is applied during the opti- 
mization of the selection. The analysis is then performed 
in a larger window 1.80 < < 1.92 GeV/c 2 . The high 
statistics decays B — > J/tp K* with the same J/ipKir fi- 
nal state are used as a control sample to evaluate the 
possible differences between data and MC. These are se- 
lected with requirements similar to those of the signal, 
except for an ttlk-k range of [0.79, 0.99] GeV/c 2 . The 
J/tp and K® candidates are constrained to their nominal 
masses [lti ] to improve the resolution of the measurement 
of the four-momentum of their parent-i? candidate. 

Candidate B mesons are formed from J/ip and D can- 
didates. Two kinematic variables are used to further re- 
move incorrectly reconstructed B candidates. The first 
is the difference AE = E* B - E£ eam between the B- 
candidate energy and the beam energy in the T(45 l ) rest 
frame. In the absence of experimental effects, correctly 
reconstructed signal candidates have AE — 0. The AE 
resolution is 7.5 MeV. For the signal region, AE is re- 
quired to be in the range [—15, +12] MeV. The second 
variable is the energy-substituted mass toes = (E^ am — 
plj 2 ) 1 / 2 , where p* B is the momentum of the B candidate in 
the T(4S) rest frame. The energy substituted mass Toes 
peaks at the nominal B mass of 5.279 GeV/c 2 for the sig- 
nal. Its typical resolution is 2.5 MeV/c 2 . A requirement 
of 5.274 < toes < 5.284 GeV/ c 2 was obtained in the op- 
timization of the signal selection. The analysis is then 
performed in the window 5.2 < toes < 5.3 GeV/ c 2 . If 
more than one B candidate is found in an event, the one 
having the smallest |A2S| is retained. 

Non-D B — ► J/ipKir decays that have toes, AE, and 
mg-vg- distributions similar to those of the signal are 
found to be the dominant contribution to the remaining 
background after selection cuts are applied. Signal events 
can be separated from non-D events by their peaking at 
the D invariant mass in the mx-n spectrum. In MC sam- 
ples, this spectrum shows a small but significant number 
of true D mesons: a D meson from the decay of one B was 
combined with a J/tp meson from the decay of the other 
B. We subtract this combinatorial background using the 
toes distribution: the rriKix distribution of the events in 
the sideband (5.21 < toes < 5.27GeV/c 2 ) is subtracted 
from the distribution of the events in the signal region, 
with a scaling factor R that is the ratio of the combinato- 
rial background in the signal region and in the sideband. 
The value of R is obtained from the integrals of the AR- 
GUS shape ^3 m n t s °f the to es distribution with a 
Gaussian function for the signal and an ARGUS shape 
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for the combinatorial background (Fig.[5J). The mss sig- 
nal window for the data is shifted by +0.6 MeV/c 2 with 
respect to MC after such a shift is observed on the J/tp 
K* control sample. We obtain R = 0.093 ± 0.011 and 
R = 0.131 ± 0.038 for J/ipD° and J/ipD+, respectively. 




m PS (GeV/c ) 



No significant signal for B — > J/tpD is observed. The 
numbers of events obtained are —0.6 ± 1.2 ± 0.2 (J/0Z3 ) 
and 1.2 ± 1.9 ±0.2 (J/0-D + ), where the first uncertainty 
is statistical, and the second one is the systematic con- 
tribution due to the uncertainties in the scaling factor 
for background subtraction, and of the D mass and mass 
resolution used in the fit. The branching fractions are: 



B 



S 



N evt x e x b 



(1) 



where 5* is the number of signal events obtained from 
the fit, N evt = 124 x 10 6 is the number of BB events in 
the data sample, and b is the product of the branching 
fractions of the secondary decays (Table HJ). 



TABLE I: Number of signal events, efficiency, secondary 
branching fraction, measured branching fraction (B) and up- 
per limit (UL) at 90% C.L. 



FIG. 2: Distribution of jties for data events with rtiKir in the 
range [1.8, 1.92] GeV/c 2 . (a): J/ipD°, (b): J/ip D+ . The fits 
are described in the text. The peaks at the nominal B mass 
are due to non-D events. 

The combinatorial-background-subtracted mx^ distri- 
butions (Fig.[3J) are fitted with the combination of a linear 
background (two free parameters) and a Gaussian signal 
in which the number of signal events S is a free parame- 
ter. The central value and the resolution of the D peak 




S e b B UL 

(%) (i(r 3 ) (io~ 5 ) (ltr 5 ) 

J/ip D u -0.6 ± 1.2 23.3 ± 0.3 4.49 -0.46 ± 0.93 1.3 

J/ip D + 1.2 ± 1.9 22.6 ± 0.2 1.15 3.7 ± 5.9 12.3 



Additional contributions to the systematic uncertainty 
of the branching fraction are described in the following. 
The relative uncertainty in the number of BB events is 
1.1%. The secondary branching fractions and their un- 
certainties are taken from PDG 16]. Other estimated 
uncertainties are from tracking efficiency (1.3% per track 
added linearly), K® reconstruction (2.5%), PID efficiency 
(3.0%) and the statistical uncertainty in the selection ef- 
ficiency. The uncertainty in the selection efficiency due 
to the uncertainty of the MC/data difference of the cen- 
tral value and of the width of the peaks in mg+g- , tueS) 
and AE is estimated from the J/4> K* control sample. A 
summary of the multiplicative contributions to the sys- 
tematics can be found in TableHTI The ratio of B° to B + 
production in T(AS) decays is assumed to be unity. The 
related uncertainty is small and is neglected here. 



FIG. 3: Background-subtracted mxi distributions, ((a): 
J/ipD°, (b): J/ipD + ) The fits are described in the text. 

are fixed in the fit. The resolution measured on the sig- 
nal MC sample is used. The agreement of data and MC 
samples has been studied with a D-meson control sample 
obtained with the same selection as for the D candidates 
of the J/tpD events. The resolutions are found to be sim- 
ilar in the data and in the MC samples, and a shift of 
— 0.6 ± 0.2MeV/c 2 of the central value is observed and 
accounted for. 



TABLE II: Summary of the contributions to the relative sys- 
tematic uncertainty (%). 







) J/^D + {Kl 7T+) 


B counting 


1.1 


1.1 


Secondary BF's 


2.7 


6.8 


Tracking 


5.2 


3.9 


K° s 




2.5 


PID 


3.0 


3.0 


MC statistics 


1.5 


1.0 


Sample selection 


1.0 


0.8 


Total 


6.9 


8.9 
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We obtain upper bounds on the branching fractions at 
90% confidence level (C.L.) assuming Gaussian statistics 
for the statistical uncertainties and taking into account 
the systematic uncertainties. We have used a Bayesian 
method with uniform prior for positive BF values in the 
derivation of these limits. We obtain upper limits of 
1.3 x 1(T 5 for B° -> J/tp D° and 1.2 x 1(T 4 for B+ 
— > J/ifi D + . The results for the neutral decay are signif- 
icantly lower than the amount needed to explain the ex- 
cess of low momentum J/i/j's after Chang and Hou flfij |. 
Together with the small upper limits on the branching 
fraction for decays B — > J/ipDir \£t\. we conclude that 
intrinsic charm as the explanation of the low momentum 
J/ip excess in B decays is not supported. 
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